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摘 要 :提出 了 一 种 基于 退 磁场 理论 的 磁 芯 探头 涡流 检测 信号 简化 数值 模拟 方法 ,基于 已 有 FEM- 
BEM 涡流 检测 数值 方法 开发 了 计算 程序 ,通过 与 商用 软件 比较 验证 了 其 有 效 性 。 首 先 ,分 析 了 磁 芯 
磁化 过 程 中 内 部 退 磁 场 的 特点 ,提出 了 磁 芯 退 磁 场 和 磁化 场 的 计算 方法 ;其 次 ,提出 了 磁 芯 线圈 激 
励磁 场 计 算 方法 ,考虑 了 激励 线圈 传导 电流 和 磁 芯 等 效 磁 流 ,利用 磁 流 环 电流 模型 和 毕 奥 萨 法 尔 定 
律 进 行 激励 磁场 计算 ;最 后 ,采用 和 激励 酸 场 计算 类 似 的 思路 进行 了 检测 信号 分 析 , 即 用 感 生 涡流 
二 次 磁场 替代 激励 线圈 磁场 来 计算 检测 线圈 磁 芯 的 退 磁场 和 磁化 场 , 并 基于 磁 芯 内 部 磁 通 变化 计 
算 检 测 电 压 信 号 。 基 于 所 提 理 论 , 本 研究 在 已 有 FEM-BEM 混合 法 涡流 检测 数值 模拟 软件 的 基础 
上 ,开发 了 可 考虑 磁 芯 影响 的 涡流 检测 数值 计算 程序 。 对 于 不 同 参 数 条 件 下 的 磁 芯 探头 涡流 检测 
问题 ,本 研究 方法 与 全 域 有 限 元 商用 软件 计算 结果 一 致 ,但 本 研究 方法 原理 更 为 直观 ,验证 了 方法 
和 程序 的 有 效 性 。 
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A numerical simulation method for ECT of probe with magnetic core 
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Abstract: A numerical tool to simulate the ECT signals of a probe with magnetic core is indispensable. 
However,the conventional full domain 3D FEM numerical method is hard to treat the probe scan and the 
numerical results cannot give clear and direct instruction for the probe design. In this paper, a numerical 
scheme taking the demagnetization field in the magnetic core into account is proposed to tackle with 
problems related to the magnetic core, and its feasibility is investigated through comparing with a 
conventional full domain FEM software. First , the features of the demagnetization field in the magnetic core 
are discussed, and methods to calculate the demagnetization field and the core magnetization. are 
developed. Second, a simulation method for the excitation magnetic field of a cored ECT probe is 


proposed. In practice, the excitation magnetic field of the ECT probe with magnetic core is obtained 
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through superposition of the magnetic fields due to the driving current in excitation coil and the equivalent 


magnetic current in core, while the magnetic current is treated with a solenoidal coil model. Third, a 


numerical scheme similar to the method for excitation coil is adopted to calculate the output signal of the 


cored pickup coil. In order to obtain the magnetization perturbation in the core due to eddy current , only 


the driving magnetic field in the magnetic core is necessary to be replaced by the field due to the eddy 


current. The pickup signal, then, can be obtained from the magnetic flux variation in pickup coil based on 


the Faraday's law. Based on the ideas above ,the conventional FEM-BEM hybrid code of A- formulation is 


updated and the ECT signals due to probes of different magnetic core configurations are calculated. The 


numerical results are found in good consistence with those of a conventional software of full FEM, which 


demonstrates the validity of the proposed method and the corresponding numerical codes. 


Key words: eddy current testing; magnetic core; numerical simulation; FEM-BEM hybrid code; 


demagnetization field 


磁 芯 是 一 种 具有 高 磁 导 率 、 低 电导 率 的 软 磁 材 
料 , 用 于 限制 和 引导 磁 通 ,广泛 用 于 变压器 .电机 、 传 
感 器 .磁头 等 机 电 设备 。 磁 芯 主 要 由 软 铁 磁 金属 
(如 纯 铁 ) 或 铁 磁性 化 合 物 (如 铁 氧 体 ) 制 成 , 因 高 磁 
导 率 特性 磁力 线 易 集 中 于 世 材 , 磁 世 中 的 磁场 强度 
会 得 到 显著 增强 , 且 由 于 软 磁 特 性 对 交 变 磁场 也 不 
会 产生 磁 滞 损 耗 。 磁 芯 线圈 的 磁场 由 载 流 线 圈 产 
生 , 虽 交 变 激励 电流 可 能 在 磁 芯 中 感 生 涡 流 , 但 由 于 
磁 芯 材料 电导 率 通常 很 小 ,其 中 的 感应 涡流 通常 可 
以 忽略 。 

涡流 检测 (eddy current testing, ECT) 是 一 种 基 
于 电磁 感应 原理 检测 金属 部 件 表 面 和 近 表 面 缺 陷 的 
无 损 检测 方法 ,是 五 大 常规 无 损 检 测 方法 之 一 ,具有 
广泛 工业 应 用 ” 。 为 增强 涡流 检测 探头 的 灵敏 度 ， 
笛 常 在 线圈 中 央 设 置 磁 忆 来 增强 激励 线圈 的 磁场 和 
检测 线圈 的 信号 。 直 观 上 设置 磁 芯 后 线圈 的 自 感 放 
大 倍数 应 接近 磁 芯 材料 相对 磁 导 率 的 大 小 ,如 铁 氧 
体 磁 芯 应 可 增加 近 千 倍 。 研 究 发 现 ,尽管 磁 世 会 增 
加 线圈 自 感 系 数 和 相应 的 涡流 检测 信号 ,但 其 增幅 
远 小 于 相对 磁 导 率 的 数值 。 为 明确 磁 芯 强化 检测 信 
号 机 理 , 指 导 ECT 磁 芯 探头 的 设计 和 信和 号 评价 , 研 
究 ECT 磁 芯 探头 的 信号 激发 和 接收 过 程 ,建立 ECT 
磁 芯 探头 检测 信号 数值 模拟 方法 ,具有 很 强 的 理论 
和 应 用 价值 。 

由 于 需 考 虑 磁 芯 的 强 磁 特 性 , 磁 世 探头 的 涡流 
检测 数值 模拟 通常 需要 采用 全 场 有 限 元 方法 ,对 磁 
芯 ` 线 圈 .检测 对 象 和 空气 区 域 采 用 不 同 的 材料 分 
区 ,其 边界 条 件 处 理 比 较 复杂 。 同 时 , 由 于 全 域 
有 限 元 方法 需要 对 线圈 和 磁 忆 进行 有 限 单 元 人 刨 分 ， 
处 理 探头 扫 查 问题 时 每 个 扫描 点 需要 重新 划分 网 


格 ,导致 计算 过 程 复 杂 和 计算 量 庞大 。XIE 等 提 
出 了 一 种 基于 退化 磁 矢 位 的 数值 模拟 方法 ,将 整个 
计算 模型 分 为 线圈 、 磁 芯 .检测 对 象 3 个 区 域 。 首 先 
求 取 基于 空心 线圈 产生 的 初始 涡流 场 ,其 次 基于 涡 
流 场 与 激励 线圈 电流 的 总 磁场 计算 磁 芯 磁场 ,然后 
修正 激励 场 再 次 计算 试 件 涡流 场 直至 收敛 。 这 种 方 
法 可 以 获得 涡流 检测 磁 芯 探头 激发 的 涡流 场 和 相应 
检测 信号 ,但 其 迁 代 计 算 过 程 相对 复杂 且 计 算 资 源 
耗费 大 ,特别 是 不 能 直观 给 出 磁 芯 对 涡流 场 和 检测 
信号 的 影响 规律 ,不 便 有 效 支撑 涡流 检测 探头 的 构 
型 设计 和 优化 。 基 于 4- 公式 的 涡流 场 有 限 元 - 边 
界 元 ( FEM-BEM ) 混合 计算 方法 和 程序 中 可 有 效 
处 理 空心 线圈 和 非 磁性 工件 的 涡流 检测 问题 ,但 由 
于 边界 条 件 问题 无 法 用 于 磁性 材料 ECT 模拟 。 本 
研究 基于 强 磁性 材料 磁化 的 磁 荷 / 磁 流 理论 模型 , 提 
出 了 一 种 新 型 简便 的 磁 世 影响 分 析 处 理 方法 ,并 将 
其 应 用 于 现 有 FEM-BEM 涡流 检测 数值 模拟 程序 ， 
有 效 实现 了 ECT 磁 芯 探头 的 激励 磁场 和 检测 信和 号 
的 分 析 计 算 。 通 过 与 传统 全 域 有 限 元 方法 结果 进行 
比较 ,验证 了 所 提 方 法 和 相应 程序 的 有 效 性 ,为 实际 
磁 芯 探头 设计 和 优化 及 信和 号 反 演 黄 定 了 基础 。 


1 理论 模型 和 方法 


1.1 表征 材料 磁化 特性 的 磁 流 模型 和 磁 荷 模型 


描述 材料 磁化 的 物理 模型 主要 分 为 磁 流 模型 和 
磁 荷 模型 “” 。 磁 性 材料 中 的 电子 自转 由 于 量子 效 
应 可 能 出 现 相 同 取向 ,导致 局 域 电 子 群 总 体 对 外 表 
现 出 磁性 。 如 将 这 些 相 同 自 旋 取 向 的 电子 看 成 一 个 
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整体 ,从 其 对 外 产生 磁场 的 观点 可 将 其 表征 为 环形 
电流 , 即 分 子 电流 。 在 没有 外 加 磁场 时 ,磁性 体 中 的 
分 子 电流 取向 随机 ,对 外 宏观 不 显 磁性 。 但 知 施 加 
外 磁场 , 磁 矩 作用 会 使 分 子 电流 趋向 同一 方向 ,导致 
对 外 宏观 呈现 磁性 , 即 在 磁性 物体 中 产生 磁化 M. 
由 于 是 利用 分 子 电流 来 表征 磁化 ,上 述 称 为 磁 流 模 
型 。 磁 流 和 物理 意义 上 的 传导 电流 不 同 , 是 为 了 表 
征 磁化 对 外 产生 磁场 特性 参照 毕 奥 萨 法 尔 定律 导入 
的 虚拟 量 。 磁 性 材料 内 部 分 布 的 分 子 电流 相互 车 加 
的 结果 会 形成 体 磁 流 , 体 磁 流 密 度 和 磁化 强度 的 关 
JJ, — Vx M 。 因 此 ,如 果 磁 忌 均 匀 磁 化 ,其 内 
部 的 体 磁 流 为 零 。 但 即使 磁 世 磁化 均匀 ,在 磁性 体 
表面 由 于 几何 不 连续 也 会 有 磁 流 的 存在 , 即 面 磁 流 ， 
其 强度 向 量 为 K, = M x n, rh n Jy loss e ri 
外 法 向 量 。 

磁体 模型 是 对 比 电荷 产生 电场 的 现象 而 导入 的 
表征 磁场 来 源 的 磁化 模型 ,其 核心 是 认为 强 磁 材 料 
中 的 磁化 对 应 于 由 正 负 磁 荷 构成 的 磁 偶 极 子 。 在 介 
质 未 磁化 时 ,各 磁 偶 极 子 取 向 杂乱 无 昔 , 其 磁 偶 极 矩 
相互 抵消 ,对 外 宏观 无 磁性 。 但 当 施加 外 部 磁场 后 ， 
由 于 各 偶 极 子 受 磁 矩 作用 会 转向 磁场 方向 ,导致 对 
外 宏观 表现 出 磁性 。 磁 性 材料 内 部 的 体 磁 荷 也 与 磁 
化 的 分 布 相关 ,其 体 密度 p。= V* M , 即 对 于 均匀 磁 
化 的 磁 世 其 内 部 的 体 磁 荷 密度 为 零 。 但 同样 ,在 磁 
必 界 面 会 表现 出 具有 面 密度 rc，= M «n 的 分 布 面 磁 
位 。 对 于 沿 轴 向 均匀 磁化 的 柱 型 磁 芯 , 仪 在 其 两 个 
端面 表现 出 集中 的 正 、 人 负面 磁 答 。 

磁 流 模型 和 磁 荷 模型 在 表征 磁性 体 对 外 的 磁性 
作用 时 各 有 特点 , 磁 流 模型 更 适合 表征 其 产生 的 磁 
感应 强度 B( 无 源 有 旋 场 ) ,而 磁 丛 模型 则 适合 表征 
其 产生 的 磁场 强度 五 ( 有 源 无 旋 场 )。 


1.2 Hi 


部 退 磁场 


磁 芯 内 部 退 磁 场 为 其 端面 磁 荷 产生 的 磁场 。 如 
图 1(a) 所 示 , 当 棒状 磁 世 在 均匀 轴 向 磁场 H, 中 被 
磁化 后 ,其 两 端 会 呈现 N.S 磁极 特性 ,携带 正 负面 
fp, WD 1 所 示 ,两 端面 磁 荷 在 磁 忆 介质 内 部 会 
产生 一 个 与 外 加 磁场 相反 的 附加 磁场 豆 " ,会 削弱 外 
加 磁场 , 称 为 退 磁场 。 为 此 磁 世 内 部 各 点 的 总 磁场 
强度 五 为 外 加 磁场 H, 和 附加 场 HH' 的 矢量 着 加 , 即 


H(r) =H,(r) +H (r) (1) 
磁 蕊 内 部 任意 点 的 磁感应 强度 为 


B=u (M+H)= m(x +1)H 
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其 中 x 为 磁化 系数 且 x +1 = 从 即 相 对 磁 导 率 。 
虽然 如 图 1(b) 所 示 磁 芯 合 成 磁场 强度 较 直 接 
施加 的 外 部 磁场 强度 小 很 多 ,但 强 磁性 磁 忆 材料 的 X 
具有 很 大 数值 , 磁 世 中 的 磁感应 强度 如 仍 较 无 磁 世 
时 大 ,可 使 线圈 内 总 磁 通 量 增 加 ,起 到 强化 线圈 检测 
信号 的 作用 。 由 于 涡流 探头 激励 磁场 为 激励 线圈 电 
流 和 磁 忆 磁 流 共同 作用 的 结果 ,同样 较 无 磁 必 时 会 
有 所 增强 。 
t- 


(a) 退 磁场 示意 图 
图 1 磁 世 内 部 磁场 强度 空间 分 布 


Fig.l Spatial distribution of the magnetic field 


(b) 磁场 强度 分 布 云图 


intensity inside around a magnetic core 
部 磁化 场 的 计算 方法 


1.3.1 磁 芯 对 于 线圈 激励 磁场 的 影响 
对 于 含 磁 忆 的 激励 线圈 ,其 产生 的 激励 磁场 磁 
感应 强度 为 线圈 电流 激发 的 磁场 B a E Bit RAE T 
产生 磁场 B... AIE, BI 
B,(r) =B oar) * Bu(r) (2) 
式 中 ,任意 场 点 7 B5 BRI BU n HESS BS EAE 
律 进行 计算 , 即 


1.3 WU 


1 


Bar) = N Jr) x V' Ir (3) 
-— Ho , y , 
B.,. (7) x i us JaC’) x V |r z t 
D ! 1 " 
[ick] oc 


式 中 :CO 分 别 为 线圈 区 域 和 磁 芯 区 域 ;$.。 为 
磁 芯 外 表面 ;J 和 K, 分 别 表示 磁 必 磁化 后 相应 的 体 
磁 流 和 面 磁 流 , 与 磁化 强度 分 布 关 系 如 1.1 节 所 示 。 
因此 ,要 得 到 磁 世 磁化 所 产生 的 激励 磁场 必须 首先 
得 到 磁 世 内 部 磁化 强度 M 的 分 布 。 
1.3.2 磁 芯 轴 向 均匀 磁化 时 的 退 磁场 和 退 磁 系 数 

首先 考虑 轴 向 均匀 磁化 磁 芯 的 磁场 特性 。 如 图 
2 所 示 , 当 柱 状 磁 芯 内 部 沿 轴 疝 均匀 磁化 时 , 磁 世 内 
部 磁化 强度 太 沿 z 回 且 处 处 相等 , 即 

M(r) = M,e, (5) 
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基于 相对 磁 导 率 的 定义 和 线性 磁 本 构 关 系 , 当 

磁 世 置 于 大 小 为 及 .的 均匀 轴 回 磁场 时 ,如 记 其 轴 疝 

退 磁 场 的 大 小 为 正 , 则 磁 芯 内 磁场 强度 和 磁化 强度 

的 轴 向 分 量 满足 以 下 关系 , 即 

(u, 19 - H.) = M, (6) 

TR d fo , E T 4635) 5] I EE E MS RE , H. 

磁化 沿 轴 向 故 在 磁 芯 侧面 也 不 存在 面 磁 荷 , 仅 在 磁 
芯 两 端面 存在 符号 相反 的 面 磁 人 荷 , 其 面 密度 为 

c, -M-:n-M, (7) 

由 于 磁 芯 内 部 的 退 磁场 由 磁 荷 产生 , 故 由 磁场 

库仑 定律 可 得 磁 芯 内 点 了 处 的 z 向 退 磁场 强度 万 为 


, Om 1 
H = — Po REN / 
s sı v Ir' - r| ow 
[vous] a) 
S. |r ex 


P, S, 、5, 分 别 为 磁 芯 上 下 端面 区 域 。 

由 式 (9) (10) 可 看 出 , 退 磁场 强度 与 磁化 强度 
之 间 具 有 线性 关系 。 由 退 磁 系 数 定 义 a =M/H', 可 
知 对 于 均匀 磁化 的 磁 芯 ,其 退 磁 系 数 a 可 表示 为 

1 1 

r | 

可 以 看 出 对 于 均匀 磁化 圆柱 形 磁 芯 ,a 是 一 个 
与 磁 芯 长 度 和 直径 相关 的 常数 。 当 磁 芯 长 度 较 直 径 
大 得 多 时 , 退 磁 系数 近似 为 0, 即 退 磁场 很 小 。 当 磁 
芯 直 径 较 其 长 度 大 的 多 时 退 磁 系数 接近 1, 即 内 部 
磁化 场 基本 为 零 。 综 合式 (6) 与 (8) ,可 求 得 磁 芯 内 
部 的 磁化 强度 ,进而 利用 式 (4) 求 出 磁 芯 对 外 产生 
的 磁感应 强度 。 


SO 


: ds" (9) 


Q 


司 2，” 磁 芯 内 部 均匀 磁化 示意 图 
Fig.2 Magnetic field and magnetization in magnetic core 
1.3.3 磁 芯 内 部 非 均匀 磁化 时 的 退 磁 场 计 算 
对 于 带 磁 世 的 盘 式 涡流 检测 探头 ,由 于 线圈 产 
生 的 磁场 不 均匀 , 且 磁 世 的 大 小 有 限 , 磁 世 内 部 的 磁 
化 强度 也 不 均 义 。 此 外 ,对 于 含 磁 世 的 检测 线圈 ,由 
于 被 检测 对 象 中 的 涡流 场 , 特 别 当 其 中 含有 裂纹 缺 
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陷 时 分 布 复杂 , 磁 世 被 涡流 二 次 磁场 磁化 时 不 均匀 
程度 更 大 。 因 此 , 当 处 理 磁 必 探 头 的 涡流 检测 数值 
模拟 问题 时 , 磁 世 内 的 磁化 必须 考虑 其 非 均 匀 性 。 
为 此 ,本 研究 提出 将 磁 芯 体积 进行 单元 离散 ,在 每 个 
小 单元 认为 磁化 均匀 ,进而 结合 线性 磁性 本 构 关 系 
与 磁 葵 磁场 的 计算 公式 ,建立 各 单元 磁化 强度 的 线 
性 代数 方程 组 来 求 取 磁 世 内 部 的 磁化 场 。 

对 于 涡流 探头 的 激励 线圈 ,通常 采用 盘 式 线圈 
和 柱状 磁 忆 ,由 于 轴 对 称 特性 其 磁 芯 内 部 的 磁场 主 
要 为 轴 向 分 量 。 对 于 检测 线圈 ,虽然 涡流 场 产 生 的 
二 次 磁场 非常 复杂 ,但 在 线圈 中 产生 感 生 电 压 的 为 
轴 向 磁场 分 量 ,考虑 到 问题 的 线性 可 县 加 性 ,同样 可 
以 局 限于 考虑 轴 向 分 量 。 因 此 ,本 研究 仅 考 虑 轴 向 
人 磁化 情形 , 即 假设 磁场 和 磁化 均 仅 有 轴 向 分 量 。 考 
虑 图 3 所 示 磁 芯 内 部 的 某 一 中 心 点 位 置 向 量 为 普 的 
小 单元 体 , 其 中 心 点 处 的 磁化 可 表示 为 


x 


图 3 磁 世 内 部 代表 磁 芯 小 单元 示意 图 


Fig.3 A concept diagram of cell division inside the core 

(u, - DUIS Gr) - H(r'))  M(r') (10) 

式 (10) 中 ,磁场 强度 和 磁化 强度 均 为 其 轴 向 分 

E, HART FER 以 下 同 )。 由 磁场 库仑 定律 ,可 

得 磁 芯 单元 的 磁 荷 在 任意 场 点 了 处 所 产生 的 = 向 退 
磁场 强度 为 

H'(r) = 2 "MUS e 

其 中 ,$S 为 磁 世 小 单元 体 上 下 端面 区 域 。 式 (11) 中 ， 

因 磁 芯 单 元 体 很 小 ,可 认为 其 内 部 的 磁化 均匀 , 故 忽 

略 了 式 (8) 中 的 体 磁 荷 项 上 且 面 磁 荷 也 仅 局 限于 单元 

体 的 上 下 端面 。 如 记 小 单元 体 上 下 端面 面积 为 AS， 

考虑 到 其 端面 磁 荷 密度 ec 等 于 轴 向 磁化 强度 M , 式 

(11) 又 可 表达 为 


IT， ,， 1 uv, 1 . 
H'(r) e -MGr)as(v Iri -r | y Pa k 


.dS (11) 


(12) 
其 中 :ri 和 7 分 别 为 小 单元 体 上 下 端面 中 心 点 位 置 
向 量 ;e: 为 z 向 坐标 向 量 。 进 一 步 利用 
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1 z-z 
Vse, =A 13 
r-r] z Ir -rb ( ) 


可 得 任意 小 单元 体 在 场 点 7 处 产生 的 退 磁 场 为 
Ir -rl de AT 
(14) 

对 于 激励 线圈 中 的 传导 电流 在 任意 点 + 产生 的 
磁化 磁场 , 则 可 利用 毕 奥 沙 法 尔 定律 利用 式 (3) 和 
H = BW/ 进行 计算 。 

如 果 磁 芯 被 离散 为 N 个 小 柱 单元 , 则 第 i 个 小 
单元 中 心 点 7 处 的 退 磁 场 为 所 有 小 单元 端面 磁 集 产 
生 退 磁场 的 和 , 即 


H(r) = a MCOAS[ 


N 
H (r,) = Y cM, (15) 
j=l 
e = HY Zn d T zn | (16) 
i 4m Ir, -rl Ira - ril 


基于 磁化 本 构 关 系 式 (6) ,可 得 各 磁 世 小 单元 
的 磁化 应 满足 以 下 方程 


N 


XU 1 = Hy,(i = 1,2,3,.…,N) 
j=1 Kr m 
(17) 
其 中 
Hy) = 去 | Jo (r!) x V! 1———À edv’ 
° ] 4r Ooil 0 |r' -r,| z 
(18) 


综合 式 (17) ~ (18) ,可 得 以 下 关于 各 单元 磁化 
强度 的 线性 方程 组 
Cy. My = Hoy (19) 
其 中 
€; ,l zJ 
(E j o (20) 
Ci Ted 一 D" =j 
3& (19) HRA N 个 未 知 数 的 满 秩 线性 代数 方 
程 组 ,可 采用 高 斯 消去 法 等 常规 方法 进行 求解 。 值 
得 指出 的 是 ,虽然 式 (11 ) 的 被 积分 函数 具有 场 点 和 
源 点 重合 的 奇异 点 ,但 由 于 式 (15) 中 的 源 点 位 于 各 
小 单元 端面 中 点 ,而 场 点 在 各 单元 的 中 心 不 会 重合 ， 
因此 已 有 效 避 免 了 奇异 点 问题 。 
1.4 磁 芯 线圈 激励 磁场 的 计算 
由 于 磁 荷 模型 适合 磁场 强度 的 计算 ,而 磁 流 模 
型 更 适合 磁感应 强度 的 计算 , 当 获 得 磁 世 内 的 磁化 
分 布 后 ,可 基于 磁 流 模型 利用 式 (2) 县 加 激励 线圈 
磁场 和 磁 必 磁 化 磁场 计算 带 磁 忆 线 圈 导 致 的 激励 磁 
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感应 强度 。 为 计算 磁 世 磁化 产生 的 磁感应 强度 Bu 
需 从 获得 的 磁化 强度 M 计算 磁 世 的 等 效 体 磁 流 与 
面 磁 流 。 实 际 上 由 于 每 个 小 磁 世 单元 体内 认为 磁化 
均匀 ,其 内 部 不 存在 体 磁 流 , 仅 在 绕 其 柱 面 存在 面 磁 
流 。 即 磁 世 整体 的 体 磁 流 表 现 为 磁 世 内 部 各 小 单元 
体 表面 的 面 磁 流 , 而 磁 忆 整体 的 面 磁 流 则 表现 为 磁 
世 最 外 侧 单元 和 空气 区 域 相 邻 表 面 上 的 面 磁 流 。 对 
各 个 小 单元 面 磁 流 产生 的 磁场 七 加 , 即 可 获得 磁 芒 
产生 的 总 磁场 为 

B.) = PY M [n x v' ra" 
(21) 
其 中 :5; 为 第 i 个 小 磁 芯 单元 的 侧 表面 区 域 ;7(7') 为 
小 磁 世 单元 侧 表面 ( 柱 面 ) 环 向 切 向 量 。 由 于 如 图 3 
所 示 各 小 单元 的 断面 形状 较为 复杂 ,对 各 小 单元 表 
面积 分 并 不 容易 。 考 虑 到 通常 涡流 检测 磁 忆 探头 的 
磁 世 为 圆柱 体 且 激励 线圈 为 盘 式 圆 线 轿 ,其 磁化 场 
具有 轴 对 称 特性 。 基 于 这 一 特性 ,本 研究 提出 以 下 
考虑 轴 对 称 特性 的 激励 磁场 计算 方法 。 
如 图 4 所 示 , 在 磁 世 磁化 强度 分 布 轴 对 称 的 情 
况 下 ,将 磁 芯 轴 向 、 径 向 、 环 向 分 别 均匀 划分 为 普 、 
mr 个 小 单元 。 这 时 轴 向 和 径 向 位 置 相同 但 环 向 
位 置 不 同 的 小 单元 具有 完全 相同 的 磁化 强度 。 
此 ,相对 于 1.3 节 给 出 的 完全 三 维 处 理 方法 这 时 的 
待 求 未 知 磁化 强度 数量 可 以 大 大 降低 。 
实际 上 , 当 考 虑 轴 对 称 特性 时 各 磁 芯 小 单元 的 
磁化 强度 向 量 M 可 用 0 =0 断面 各 小 磁 世 单元 的 磁 
化 强度 向 量 M RS AE o 


M - TM, 


n.Xn,Xng 


(22) 


n.Xn, 


图 4 磁 忆 场 计 算 轴 对 称 单元 示意 图 
Fig.4 Concept of axial cell division of magnetic 


core for magnetic field calculation 
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其 中 ,了 是 基于 轴 对 称 特性 的 三 维 磁化 强度 和 二 维 
轴 对 称 磁化 强度 变量 的 转换 矩阵 。 

基于 式 (19) 和 (22 ) 可 得 关于 M, 的 线性 代数 方 
程 组 为 


T'CTM, = T'H, (23) 

求解 方程 (23 ) 获得 M, 后 ,利用 式 (22) 即 可 获 
得 所 有 小 单元 的 磁化 强度 。 

在 轴 对 称 情况 下 ,存在 比 式 (21 ) 更 方便 的 方法 
计算 磁 芯 磁化 所 产生 的 激励 磁场 。 实 际 上 ,如 图 4 
所 示 , 由 于 轴 对 称 时 磁 流 仅 为 各 圆 环 单元 侧面 上 的 
面 磁 流 ,各 单元 内 部 的 体 磁 流 为 零 ,为 此 各 环 面 上 的 
面 磁 流 可 用 0 =0 断面 单元 的 磁化 唯一 表征 , 即 

! M? — MË j x n, 
Kus. 
M'!,j-n, 
2 ons PT 
AUS n, ux n RAE S 
各 小 单元 外 侧 表面 的 中 心 位 置 , 其 电流 大 小 为 
s. (25) 
HPI i RA NEN S EUM 
小 ;hj 为 相应 单元 的 轴 向 厚度 。 这 些 环 电 流产 生 的 
磁场 ,可 基于 毕 奥 萨 法 尔 定律 利用 椭圆 积分 解析 计 
算 。 实 际 上 ,如 记 任 意 场 点 P 的 位 置 柱 坐 标 为 (p,， 
0,,5,), 令 Ro=[(py+p,) * (x, -2,) ] ^ WPA 
B5 24 I Jy v. BE 46 qn] RH dl 6] y Eit I 9] 2 CX I8] 4 
量 为 零 ) 
E aS »n EX [2 Dj *p» 


(pj -pp) 十 (zp -u) 


(24) 


i + (zp -a) 


a - (zp — 
“pp (p; FSE + (zp =z) 


-Hy 24 


(27) 
其 中 :K 和 分 别 表示 第 一 类 和 第 二 类 椭圆 积分 ; 
p; 3 分 别 为 轴 向 第 i 个 、 径 向 第 j 个 环 电流 的 径 向 和 
轴 向 坐标 。 

将 式 (26) (27) 所 得 磁 芯 磁化 产生 的 磁场 和 式 
(3) 的 线圈 流 激发 的 磁场 相 加 , 即 可 得 到 带 磁 芯 激 
励 线圈 产生 的 总 磁感应 强度 , 即 激励 磁场 。 


1.5 磁 芯 探头 检测 线圈 检 出 信号 的 计算 
根据 法 拉 第 电磁 感应 定律 ,涡流 检测 的 检 出 信 
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号 是 和 检测 线圈 交 链 磁 通 的 变化 所 致 。 为 此 计算 检 
" 信和 号 需要 计算 由 于 涡流 场 在 检测 线圈 磁 芯 中 导致 
的 总 磁 通 量 , 这 时 同样 需要 首先 计算 出 由 于 涡流 场 
导致 的 检测 线圈 磁 忆 内 的 磁化 分 布 。 考 虑 到 问题 的 
类 似 性 ,这 时 只 需 将 前 述 考虑 激励 线圈 磁化 分 布 时 
的 磁化 场 现代 之 以 检测 对 象 中 涡流 场 产生 的 次 生 
磁场 了 ,i, 即 可 。 由 于 涡流 检测 对 象 特 别 是 包含 裂 
纹 等 缺陷 时 为 非 轴 对 称 , 检 测 线圈 磁 世 的 磁化 计算 
需要 采用 1.3.3 给 出 的 三 维 磁 化 分 布 计算 方法 。 

实际 上 ,根据 毕 奥 萨 伐 尔 定 律 ,检测 线圈 磁 世 
任意 点 由 于 涡流 的 次 生 磁 场 强度 为 
1 


H ar) 三 zl Jay r’) x rm 


- r| 


dv’ 


(28) 
E RA SUI 2 PEAH AE, i h 2R FET E OBS R A A A 
化 强度 计算 时 同样 可 以 只 考虑 轴 向 分 量 。 利 用 式 
(18) 〈22) 〈23 ) ,可 得 检测 线圈 磁 芯 中 由 于 涡流 场 
导致 的 磁化 密度 的 控制 方程 组 为 
CM aay = Ha, (29) 
其 中 ,五 ,3 为 涡流 场 在 检测 线圈 磁 芯 各 小 单元 中 心 
点 产生 的 轴 癌 磁场 强度 向 量 。 求 解 方程 (29) ,可 得 
检测 线圈 磁 世 中 的 磁化 强度 分 布 。 
根据 相对 磁 导 率 的 定义 和 线性 磁性 本 构 关系 ， 
检测 线圈 磁 芯 内 部 各 点 轴 向 磁感应 强度 可 由 相应 点 
的 磁化 强度 直接 计算 获得 , 即 


BG) o PM as Qr) (30) 


1 
当 激励 电流 为 角 频 率 为 w 的 正 旋 信号 时 ,由 法 
拉 第 电磁 感应 定律 , 磁 忆 内 的 磁 通 变化 导致 的 检测 
线圈 信号 为 


- Diol ^ 

inisedd eio sap nasus 
的 磁 芯 区 域 。 
图 5 为 检测 线圈 和 磁 芯 在 涡流 场 作用 下 的 磁 通 
分 布 示意 图 ,由 于 磁 芯 的 聚 磁 效应 ,和 检测 线圈 交 链 
的 磁 通 主要 为 聚集 在 磁 芯 内 部 的 二 次 磁场 。 而 在 磁 
芯 外 部 由 于 磁 芯 磁化 导致 的 磁 通 和 涡流 导致 的 磁 通 
方向 相反 ,基本 可 以 忽略 。 为 此 ,检测 线圈 总 输出 电 
压 信号 可 用 下 式 计 算 。 


n Li 
AV = jo 2 2 AS; PE Meas (32) 


M a, (r)dS (31) 
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其 中 :ASx 为 第 下 检测 线圈 位 置 的 第 堪 的 小 柱 单元 
的 端面 积 ;L/L 为 第 i 下 线圈 面 内 各 小 磁 必 单元 的 总 
数 。 对 于 绝对 式 盘 式 涡 流 检测 线圈 ,如 其 单 古 激励 
电流 为 幅 值 为 ,, 则 线圈 的 阻抗 变化 信号 为 

AQ = AV/I, (33) 


图 5 含 磁 世 检 出 线圈 磁 通 分 布 示意 图 


Fig.5 Magnetic flux distribution of cored detection coil 


2 程序 开发 和 数值 计算 结果 


2.1 磁 芯 探头 涡流 检测 FEM-BEM 计算 程序 开发 


基于 上 节 所 述 理论 和 作者 等 开发 的 FEM-BEM 
混合 法 涡流 检测 数值 模拟 软件 ,编写 了 带 磁 芯 涡 
流 探 头 检测 数值 模拟 程序 。 本 程序 无 需 对 空气 区 
域 、 磁 芯 和 线圈 进行 有 限 元 网 格 刨 分 ,不 仅 可 减 小 计 
算 区 域 , 还 能 方便 地 模拟 探头 扫 查 检测 信号 。 本 研 
究 开发 程序 模块 包括 以 下 两 个 部 分 。 

1) 磁 忆 激励 线圈 在 检测 对 象 中 的 总 激励 磁场 计 
算 。 基 于 式 (20) ,首先 开发 了 人 处理 激励 线圈 磁 芒 全 
化 分 布 的 系数 矩阵 C 的 计算 子 程序 ;其 次 ,基于 自由 
空间 激励 线圈 磁场 式 (18 ) ,开发 了 磁 必 内 磁场 强度 
向 量 有 ,的 计算 子 程序 ;最 后 ,开发 了 求解 式 (19) 获 
取 磁 世 磁 化 分 布 , 及 基于 式 (24) 获 取 环 电流 大 小 和 
位 置 的 程序 。 

2) 磁 心 检测 线圈 的 检 出 信号 计算 。 首 先 基于 式 
(26) ~ (27) 开 发 了 计算 涡流 场 在 磁 芯 小 单元 中 点 
产生 的 磁场 强度 ,形成 方程 式 (30 ) 右 端 项 吾 必 的 子 
程序 ;其 次 ,开发 了 基于 式 (15) (19) (21) JE UIS 
测 线圈 磁 芯 的 系数 矩阵 C, 并 求解 方程 (30 ) 获取 
Mu 的 子 程序 ;最 后 ,利用 线性 磁性 本 构 关 系 , 基 于 
式 (30) ~ (33) 计算 检测 信号 的 磁 芯 线圈 检测 信和 号 
计算 程序 。 
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上 述 程序 模块 利用 FORTRAN 程序 语言 编写 ， 
集成 到 了 原 有 FEM-BEM 混合 法 程序 ,实现 了 磁 世 
涡流 探头 的 检测 信号 的 有 效 计 算 。 


2.2 全 域 有 限 元 磁 芯 探头 涡流 检测 特性 分 析 


为 明确 磁 芯 特性 参数 对 磁 芯 探头 涡流 检测 信和 号 
的 影响 ,首先 利用 商用 软件 对 典型 磁 世 线圈 探头 涡 
流 检 测 问题 进行 了 全 域 有 限 元 数值 模拟 计算 。 涡 流 
检测 对 象 设 定 为 Inconel 旬 基 合金 钢 平 板 , 其 尺寸 和 
材料 参数 见 表 1。 探 头 激励 线圈 采用 圆柱 磁 芯 盘 式 
线圈 ,其 线圈 、 磁 芯 斥 才 、 检 测 条 件 和 材料 参数 如 表 
2 所 示 。 

表 1 平板 模型 检测 对 象 参数 

Tab.1 The parameters of the plate model that detect the object 
长 度 /mm 宽度 /mm 厚度 /mm 电导 率 /MSm-! 相对 磁 导 率 
40.0 40.0 1.0 1.0 1.0 

R2 线圈 参数 .检测 条 件 及 圆柱 磁 芯 参数 


Tab.2 Coil parameters detection conditions and 


cylindrical core parameters 


项 目 参数 数值 
线圈 参数 内 半径 /mm 0.6 
外 半径 /mm 1.6 
轴 向 厚度 /mm 0.8 
ri c 140 
检测 条 件 总 激励 电流 /A 0.1 
激励 频率 /kHz 300 
提 离 距离 /mm 0.5 
磁 芯 参数 直径 /mm 1.2 
高 度 /mm 0.8 
相对 磁 导 率 2 000 
电导 率 0 


全 域 有 限 元 计算 采用 四 面体 单元 ,对 磁 芯 、 线 
圈 空气 区 域 和 检测 对 象 平板 均 进 行 了 有 限 元 网 格 
刨 分 ,对 有 无 磁 芯 条 件 下 的 激励 线圈 磁场 .涡流 场 和 
检测 线圈 阻抗 信号 分 别 进 行 了 计算 。 

图 6(a) 给 出 了 磁 世 激励 线圈 产生 的 激励 磁场 
的 矢量 分 布 和 在 检测 平板 表面 的 磁场 大 小 分 布 云图 
的 全 域 有 限 元 计算 结果 。 设 置 磁 世 后 最 大 激励 磁场 
会 增加 1.8 倍 , 且 更 集中 于 线圈 下 方 。 图 6(b) 给 出 
了 绝对 式 探头 有 无 磁 世 时 计算 所 得 检测 信号 和 激励 
频率 的 关系 。 可 以 发 现 ,对 于 盘 式 线圈 绝对 涡流 探 
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头 , 加 入 磁 芯 可 使 检 出 信号 增 大 一 倍 左右 。 


(a) HANA zs R 


^ 无 磁 芯 
30F x 有 了 磁 世 


25 


2.0 


阻抗 信号 幅 值 /Q 


0 50 — 100 150 200 250 300 
激励 频率 /kHz 
(b) 阻抗 信号 与 激励 频率 关系 
6 ”商用 软件 典型 仿真 计算 结 


Fig.6 Typical simulation results of a commercial software 


Ea 


2.3 本 研究 方法 和 程序 的 有 效 性 验证 


为 验证 本 研究 提出 的 磁 芯 探头 涡流 检测 数值 模 
拟 方法 有 效 性 ,采用 本 研究 程序 对 上 述 线圈 和 检测 
对 象 的 涡流 检测 信号 进行 了 计算 ,并 与 全 域 有 限 元 
计算 结果 进行 了 比较 。 计 算 中 检测 对 象 平板 在 x 方 
向 划分 为 17 个 单元 y 方向 划分 为 28 个 单元 .z 方 
向 (厚度 方向 ) 划 分 为 4 个 单元 (合计 1904 个 六 面 
体 单元 和 1 132 个 矩形 边界 元 ) 。 为 考察 探头 对 缺 
陷 的 检测 灵敏 度 , 在 平板 的 中 央 设 置 了 沿 y 向 的 长 
10 mm 、 宽 0.2 mm .深度 为 板 厚 50% 的 表面 裂纹 , 计 
算 了 探头 在 x = 0 、 提 离 0. 5 mm 位 置 沿 方向 扫 查 
20 mm 的 检测 信号 ( 扫 查 间隔 1 mm) 。 涡 流 探头 选 
a ea a 

磁 芯 线圈 。 计 算 模型 和 有 限 元 网 格 如 图 7 所 示 。 

为 验证 方法 的 有 效 性 ,本 研究 对 A B,C 共 3 种 
不 同 尺寸 磁 芯 进行 了 数值 计算 ,各 磁 芯 参数 如 表 3 
所 示 。 计 算 中 磁 世 在 轴 向 、 径 向 和 环 向 分 别 均匀 划 
4398 n, 2 10 mw =10 n, =12 个 单元 。 激 励 场 的 计算 
分 别 采用 了 三 维 模型 计算 磁化 分 布 和 轴 对 称 模型 计 
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算 激励 磁场 强度 。 检 测 信号 的 计算 则 采用 三 维 模型 
计算 磁化 ,利用 式 (32) 计算 检测 信号 。 


(a) 涡流 检测 缺陷 扫描 示意 图 


A, 


(b) 计算 程序 中 检测 试 件 网 格 划分 


图 7 数值 计算 模型 和 有 限 元 网 格 
Fig.7 Numerical model and finite element mesh 

全 域 有 限 元 计算 采用 了 2.2 节 所 示 计 算 条 件 ， 
并 通过 移动 裂纹 计算 裂纹 扫 查 检测 信号 。 获 得 全 部 
21 个 扫 查 点 的 涡流 检测 信号 ,本 研究 程序 和 全 域 有 
限 元 的 计算 时 间 分 别 为 0.13 h 和 1.5 4b, 本 研究 程序 
具有 更 高 的 计算 效率 。 

表 3 圆柱 磁 芯 模型 尺寸 


Tab.3 Cylindrical core model size 


模型 参数 数值 


标准 磁 芯 : 磁 芯 人 A 长 度 /mm 0.8 
直径 /mm 1.2 
长 型 磁 世 : 磁 芯 B 长 度 /mm 1.6 
直径 /mm 1.2 
细 径 磁 芯 : 磁 世 C 长 度 /mm 0.8 
直径 /mm 0.8 


1 ) 磁 心 A( 标 准 磁 世 ) 计算 结果 比较 

图 8 给 出 了 磁 世 A 激励 线圈 在 试 件 上 表面 产生 
的 激励 磁场 分 布 计算 结果 。 其 中 图 8(a) (pb) 分 别 
为 本 研究 程序 和 全 域 有 限 元 计算 磁场 大 小 的 二 维 分 
布 ,图 8(c) 为 其 在 x=0 线 上 的 大 小 比较 。 可 以 发 
现 两 者 恨 好 一 致 ,说 明了 本 研究 磁 芯 线圈 激励 场 计 
算 方法 和 程序 的 有 效 性 。 图 9 给 出 了 含 磁 世 A 线圈 
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裂纹 检测 信号 的 相 平面 计算 结果 比较 ,所 示 信 号 为 
含 裂纹 信号 计算 结果 和 无 裂纹 信号 计算 结果 的 差 
分 , 即 裂纹 导致 的 信号 扰动 。 可 以 看 出 两 者 同样 基 
本 一 致 ,进一步 确认 了 本 研究 方法 的 有 效 性 。 
—T———— D 
prr m 327.0 
292.9 
258.8 
224.7 
190.6 
156.5 
122.4 
88.3 
54.2 
20.1 
-14.0 


磁场 强度 /(A，m ，) 
— 337.0 
301.2 
2654 
229.6 
193.8 
158.0 
122.2 
86.4 
50.6 
14.8 
-21.0 


RREA * m) 


250 


磁场 强度 /(A。m ) 
A a 


[1 
e 


e 


-50 
—0.003 —0.002 


—0.001 0 
扫描 位 置 /m 


(c) 中 线 上 的 两 者 对 比 


0.001 0.002 0.003 


图 8 标准 磁 芯 线圈 激励 磁场 在 试 件 上 表面 处 分 布 
Fig.8 The distribution of excitation magnetic field 
intensity due to the standard core coil at the 


upper surface of the specimen plate 
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—0.05 0 0.05 0.10 


实 部 /Q 


图 9 标准 磁 芯 线圈 检 出 信号 对 比 


Fig.9 Comparison chart of standard magnetic 


0.15 


core coil detection signal 
2) Bit B PIRE C 的 计算 结 
图 10 ~ 11 分 别 给 出 了 采用 磁 世 B 时 的 激励 磁 
场 分 布 和 裂纹 检测 信号 的 计算 结果 ,图 12 ~ 13 分 别 
给 出 了 采用 磁 世 C 时 的 激励 磁场 分 布 和 裂纹 检测 信 
号 的 计算 结果 。 
l pc m) 


229 
203 


磁场 强度 /(A* m^) 


(a) 本 研究 结果 
磁场 强度 /(A.m ) 

257.0 
229.4 
201.8 
1742 
146.6 
119.0 
91.4 
63.8 
36.2 
8.6 
-19.0 


250] - 


200 


磁场 强度 /Am ) 
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= 
© © © 


[91 
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—0.003 -0.002 -0.001 0 0.001 0.002 0.003 
扫描 位 置 /m 
(c) 两 者 对 比 图 
图 10 磁 芯 3 线圈 激励 磁场 在 试 件 上 表面 处 的 分 布 


Fig. 10 The distribution of excitation magnetic field intensity 


of the long type core coil on the upper surface of the specimen 


0 0.05 0.10 0.15 0.20 0.25 
实 部 /9 
图 11 磁 芯 3 线圈 检 出 信号 对 比 
Fig. 11 Comparison chart oflong type magnetic 
core coil detection signal 
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图 12 RA C 线圈 激励 磁场 在 试 件 上 表面 处 磁场 分 布 比较 


Fig. 12 The excitation magnetic field intensity distribution 


0.001 0.002 0.003 


of the thin type core coil on the upper surface of the specimen 
可 以 发 现 ,本 研究 方法 和 程序 对 不 同 磁性 参数 
均 可 获得 和 全 域 有 限 元 很 好 一 致 的 计算 结果 ,显示 
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了 方法 的 有 效 性 和 适应 性 。 同 时 ,比较 图 9、11、13 
的 裂纹 检测 信号 ,不 难 发 现 更 大 的 磁 芯 长 度 和 填充 
率 对 提高 检测 灵敏 度 有 效 。 


实 部 /Q 


图 13 Cit C 线圈 检 出 信号 对 比 


Fig.13 Comparison chart of thin type core magnetic 


coil detection signal 
综 上 ,本 研究 方法 和 程序 与 全 域 有 限 元 软件 计 
算 结 果 很 好 一 致 ,利用 FEM-BEM 混合 方法 有 效 实 
现 了 含 磁 忆 涡 流 检 测 问题 的 数值 模拟 , 且 本 研究 方 
法 无 需 对 空气 区 域 进 行 网 格 划分 ,更 便于 进行 扫 查 
言 号 的 模拟 计算 ,具有 更 高 的 效率 。 


3 结 it 


基于 磁 荷 / 磁 流 模型 磁化 理论 ,本 研究 提出 了 仿 
芯 涡 流 检测 探头 的 激励 磁场 和 检测 信号 的 计算 广 
法 ,推导 了 相应 的 理论 公式 ,实现 了 基于 FEM-BEM 
混合 法 的 磁 芯 探头 涡流 检测 信号 计算 ,解决 了 全 域 
有 限 元 计算 带 来 的 计算 效率 和 直观 性 问题 。 论 文 得 
出 了 以 下 主要 结论 。 

1) 基 于 磁 荷 / 磁 流 模型 的 方法 可 有 效 实现 含 磁 
芯 激 励 线圈 的 激励 磁场 和 合 磁 芯 检 测 线圈 检测 信号 
计算 。 

2) 本 研究 方法 和 程序 的 计算 结果 与 全 域 有 限 元 
计算 结果 一致 ,由 于 无 需 对 线圈 和 磁 芯 进行 有 限 元 
网 格 蚀 分 ,本 研究 程序 具有 更 好 的 计算 效率 和 适 
应 性 。 

3) 对 于 常规 盘 式 涡流 检测 探头 , 磁 芯 会 提高 检 
测 信号 , 较 长 的 磁 芯 和 较 大 的 磁 芯 填充 率 具有 更 好 
的 裂纹 检测 性 能 。 
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